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We examine linear and quasi-liner stages of Cherenkov-drift instability in the relativistic magne¬ 
tized electron-positron plasma. The external magnetic held lines are assumed to be slightly curved. 

In this case the curvature drift of relativistic beam particles plays decisive role in the development of 
the instability. Quasi-linear relaxation of the relativistic beam leads to diffusion of the resonant par¬ 
ticles in the momenta space. The expressions for diffusion coefficients of Cherenkov-drift instability 
are obtained. 


I. INTRODUCTION 


Gherenkov-drift instability was suggested by Kazbegi, 
Machabeli & Melikidze Sill as a possible mecha¬ 
nism for generation of pulsar radio emission and later it 
was approved in Q. In those works the linear theory of 
Cherenkov-drift instability was developed. It was shown, 
that in the pulsar magnetosphere due to Cherenkov-drift 
instability the orthogonally polarized plasma waves are 
excited. These waves can escape from the magnetosphere 
and reach observer as a pulsar radio emission. The nec¬ 
essary condition for the development of Cherenkov-drift 
instability (as for an usual Cherenkov instability) is a 
presence of a beam of particles in the relativistic magne¬ 
tized pair plasma (consisted of relativistic electrons e~ 
and positrons e’*'). 

Generally, Cherenkov type instabilities develop due to 
a resonant interaction between waves and particles of a 
beam. The resonance occurs, when the electric field vec¬ 
tor E and the wave vector k of generated waves have 
got components along direction of the beam velocity v 
(E • V 7 ^ 0 and k • v 7 ^ 0). So, in the magnetized plasma 
the transverse waves (E _L B _L k) propagating along the 
external magnetic field ( Bq || k ) cannot be generated 
by an usual Cherenkov instability (because it develops 
on the beam particles moving along the external straight 
magnetic field lines: v |j Bq and E • v = 0 ). 

Cherenkov-drift instability develops when the beam 
particles move along slightly curved magnetic field 
(SCMF) lines and, hence, drift across the plane where 
the curved lines lie. The drift motion of beam parti¬ 
cles provokes generation of purely transverse as well as 
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longitudinal-transverse waves. 

Generally there are two most important effects caused 
by the the particle relativistic motion along the SCMF 
line: curvature drift and curvature radiation. The drift 
velocity is directed across the plane of the SCMF lines 
and is given by the following expression: 


UJbRb 


( 1 ) 


Here Ud denotes the drift velocity of electrons (positrons 
drift with the same velocity but in the opposite direc¬ 
tion); ujb = eBimc is the cyclotron frequency of elec¬ 
trons; Rb is the curvature radius of the magnetic field 
line; 7 is Lorentz factor of a particle; c is speed of light 
and U|| is the component of v along the magnetic field 
line. If 7 » 1, the value of drift velocity Ud could be 
signihcant. 

A single particle, moving along the curved held line, 
radiates so called curvature radiation which can be eas¬ 
ily described as a synchrotron radiation in an effective 
magnetic held (see e.g. |^). In 1975 Blandford inves¬ 
tigated the curvature radiation of plasma howing along 
the SCMF lines. The problem was studied in the limit 
of inhnite magnetic held Bq 00 and it was shown that 
there is no radiation at all: the waves, radiated by each 
particle, are absorbed by another one. This result was 
conhrmed later in papers g gl, 0, 0, where spa¬ 
tially unlimited plasma how was considered. Then in the 
paper by Asseo, Pellat and Sol 0 the sharp boundary 
was assumed at the edge of the how propagating along 
the curved held line and the possibility of the waves ex¬ 
citation was shown at this boundary. If plasma how has 
zero width the instability is reduced to that of Goldreich- 
Keeley [0. 

Development of Cherenkov type instability, taking into 
account the curvature drift motion, was studied and 
growth rate was calculated in g g g 0, 0, 0, 0, 0 
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for different particular cases. These results were con¬ 
firmed later after thorough investigation of the problem 
in &|ii- The instability was called a Cherenkov-drift 
instability. 

Presence of the curved magnetic field lines is the neces¬ 
sary condition for both the Cherenkov-drift radiation in 
plasma and the single particle curvature radiation in vac¬ 
uum. However, the Cherenkov-drift radiation still could 
not be interpreted as a plasma curvature radiation, anal¬ 
ogous to the single particle curvature radiation: as it will 
be shown below, the Cherenkov-drift radiation is not gen¬ 
erated in the case of Ud 0. On the other hand, it is evi¬ 
dent that the drift velocity equals to w 0 if we assume 
that Hq —> oo. However a single particle radiates even for 
the infinite intensity of magnetic field. Moreover, the sin¬ 
gle particle radiates the vacuum wave, while the proper 
waves of the medium (i.e. relativistic electron-positron 
plasma) are generated in the case of Cherenkov-drift in¬ 
stability. This particular point was not considered in the 
works by Blandford 0 and Melrose ||]. Polarization of 
these waves strongly differs from that of vacuum waves. 

Brief examination of the linear theory of Cherenkov- 
drift instability is discussed in section 2. In section 3 the 
quasilinear equations for the Cherenkov-drift instability 
are obtained. In section 4 coefficients describing diffusion 
of particles in momentum space are evaluated. Alteration 
of plasma distribution function is studied. The results are 
summarized in section 5. 


II. THE LINEAR THEORY 

Properties of relativistic electron-positron plasma were 
carefully investigated in series of works |2^, |2^, 
and are still the subject of concern The electron- 

positron plasma of pulsar magnetosphere consists mainly 
of two components: a bulk of plasma particles with high 
density and low Lorentz factors (qp ~ 3-^10) and a beam 
of the primary particles ejected from the stellar surface. 
Density of the beam rib ncj, where ncj is so called 
Goldreich-Julian density. Density of secondary plasma 
Tip ~ Ti-rib, where E is the Sturrock ’multiplication factor’ 
(S ~ 102 ^ 10^) @. 

The electron-positron plasma differs from the electron- 
ion plasma by lack of gyrotropy. Consequently, spectra 
of waves propagating in the e“e^ plasma is simpler than 
that in the electron-ion plasma. It consists of only four 
types of waves which correspond to four branches on the 
diagram tu(fc), where w is the wave frequency (see Fig.^. 

One of them is high-frequency transverse electromag¬ 
netic wave totally located in superluminal area. Its phase 
velocity Vph = uj/k exceeds speed of light c, hence, it is 
not of interest in our discussion below. 

The second branch represents the dispersion of purely 
transverse linearly polarized electromagnetic wave. We 
call this dispersion curve t mode. It is totally located 
in subluminal area and therefore could be generated by 


beam particles. Its electric field vector E* is perpendic¬ 
ular to the plane of wave vector and external magnetic 
field (k,Bo). 

The remaining two dispersion curves on the u;(/c) dia¬ 
gram describes the longitudinal-transverse waves propa¬ 
gating in relativistic e“e^ plasma. One of them is almost 
superluminal. This wave is purely longitudinal if it prop¬ 
agates strictly along the magnetic field line (k || E || Bq), 
and, in this case called Langmuir wave associated with 
longitudinal oscillations of the charge density. If an an¬ 
gle 'd between k and Bg increases, the component of wave 
electric field E starts to grow across k: Langmuir wave 
transforms to the longitudinal-transverse wave denoted 
in Fig.|^ as lt 2 mode. If the angle 'd is small enough 
~ where ujp = sjAirripe^jm is the 

Langmuir frequency), lt 2 mode is almost longitudinal 
and crosses uj = kc line. In this case, U 2 mode could 
be excited, if Cherenkov resonance condition oj = fc||U|| 
is fulfilled. However, for the resonant particles of pri¬ 
mary beam, growth rate of the instability is very small 
. The wave leaves from the interaction area so quickly 
that no time is left for significant growth of the wave am¬ 
plitude. In the case of oblique propagation, 'd > 'do, lt 2 
mode is totally superluminal. Therefore, lt 2 mode could 
not be generated at all by particles of beam. 

Another longitudinal-transverse wave is denoted in 
Fig.|^ as Iti mode. This mode, like t-wave, is located to¬ 
tally in subluminal area and can easily be generated by 
plasma particles. Its electric field vector E*‘ is located in 
(k,Bo) plane. Iti mode is vacuum wave, if it propagates 
along the magnetic field lines (k I Bg). Its dispersion 
curve merges with t mode (see Fig.n ) and can be arbitrar¬ 
ily polarized. In the case of oblique propagation, electric 
field of Iti wave has the component E||‘ along the exter¬ 
nal magnetic field, thereby involving plasma particles in 
longitudinal oscillations. 

Generation of /ti-mode, propagating in perpendicular 
direction to the plane of SCMF lines, is connected with 
the drift motion of the particles. These waves are also 
known as drift waves |^, ^ . 

It should be mentioned that, while describing the 
waves in relativistic e“e“'" plasma, some authors some¬ 
times use terminology which, in our opinion, appears to 
be misleading. For example, since the work by Arons and 
Barnard p^ , the superluminal longitudinal-transverse 
wave {lt 2 mode) was called an ordinary (O) mode, the 
subluminal transverse wave {t mode) - an extraordinary 
(A) mode and the subluminal longitudinal-transverse 
wave {Iti mode) - an Alfven mode. However ’ordinary’ 
and ’extraordinary’ are generally related to the waves 
propagating across the external magnetic field in usual 
electron-ion plasma Moreover, t wave (so called 

X mode) is the purely transverse wave and its analogue 
does not exist in the electron-ion plasma. In electron- 
positron plasma the terminology ordinary and extraordi¬ 
nary are used for the waves propagating almost along the 
magnetic field. As for the Alfven mode, in electron-ion 
plasma such a name is used for an almost linearly po- 
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FIG. 1: Dispersion curves of the waves in the relativistic magnetized electron-positron plasma. The solid line corresponds 
to oj = fee. It divides the plane {ui,kc) into superluminal and subluminal areas. The high-frequency branch of t waves (in 
the superluminal area) is dehned as HF mode. Low-frequency and high-frequency modes of It waves are defined as Iti and lt 2 
modes respectively. 


larized, transverse electromagnetic wave with frequency 
Lo (where tus. is the cyclotron frequency of ions). 

In the case of fc ^ oo and to this wave trans¬ 

forms into right-hand polarized ion-cyclotron mode with 
frequency u> ^ ojBi and left-hand polarized electron- 
cyclotron wave with frequencies uj ss ujb^- The latter 
case corresponds to the fast magnetosonic wave and is 
called helicon for frequencies tus. < to < ojb^- There¬ 
fore, in relativistic e“e+ plasma, we prefer to call the 
dispersion curves t, lt 2 and Iti modes respectively, hence 
avoiding a possible confusion with dispersion curves in 
electron-ion plasma. 

To obtain expression for dielectric permittivity tenor 
components ey (u;,k) in the case of SCMF lines, it is 
handy to consider the problem in cylindrical coordinates 
{x, r, ip) and direct x axis perpendicularly to the 
plane of the curved magnetic field lines (see Fig.®. In the 
papers by Kazbegi, Machabeli and Melikidze g ||, p0[ , 
it was shown that t and It waves could be generated by 
particles of the beam when the following resonance con¬ 
dition is satisfied (for the resonant values of parameters 
io = Wo; = Vo and Ud = uq): 


Aw = w - - kxUd = 0. (2) 

Mechanism of wave generation is a modification of the 
well-known beam-plasma instability. However, this in¬ 


stability differs significantly from the usual beam-plasma 
instability Q. 

The expressions for the growth rates of the different 
waves are as follows: 

pt _ , 

2 Wk ^ ^ 

for transverse (t) wave and 


2 Wk k\ 


( 4 ) 


for longitudinal-transverse {It) wave (here ujb is Lang¬ 
muir frequency of resonant particles of the beam). These 
expressions were obtained in papers & I H- As for 
growth rate of the drift wave, it was obtained in paper 

%■ 


Td 



( 5 ) 


The reason for generation of t, It and drift waves is 
presence of the beam in the relativistic pair plasma, al¬ 
though the waves cannot be excited without drift motion 
of particles. Indeed, expressions i),® and ®) are equal 
to zero if Ud = 0. All those waves (see Eqs. are purely 
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FIG. 2: The cylindrical frame of reference (a:, r, ip) and a local Cartesian frame of reference (x, r, y). x axis is directed up 
from the plane of the figure. Bq is the external magnetic field; Rb is the radius of curvature of a magnetic held line. 


or almost transverse waves. For t and Iti waves the elec¬ 
tric field vector is perpendicular to the external magnetic 
field (see Eqs. ||,|). As for the waves generated by the 
usual beam instability, both their electric field vector E 
and the wave vector k are directed along the external 
magnetic field Bq. 

In the following section we study the quasilinear equa¬ 
tions which are significantly different from those of the 
usual beam-plasma instability. 


fao{p,fit) -f/ai(r,p,t), E(r,t) = Ei(r,t) and B(r,t) = 
Bo(r, -|-Bi(r, t). Then, averaging equation (g) over 
fast oscillations and assuming {fai) = (Bi) = (Ei) = 0, 
ifao) = fao and fao ^ fai, the following equations can 
be obtained: 


dfao 

dptt 


rriaCdp \ 1 J 



= QL, 


(7a) 


III. QUASI-LINEAR EQUATIONS 


Let us study the quasilinear theory of Cherenkov-drift 
instability, using the collisionless kinetic equation in the 
following form:, 


dfa 

dt ^ dr dp 


Qa 

rriaC 



p X B 



0 . ( 6 ) 


Here fa = /a(r, p,t) is the distribution function of par¬ 
ticles of sort a, E(r,t) and B(r,t) are the electric and 
magnetic field vectors, respectively. 

According to standard scheme, in order to obtain a 
system of quasilinear equations, distribution function, as 
well as electric and magnetic field vectors, have to be 
divided into the main and oscillating parts: fa{r,p,t) = 


dfal 

dt 


CP dfal Qa / P X Bo \ dfal 

7 dr ruaC V 7 J dp 

gg -h P ^ ^/go 

niaC V ^ 1 J dp ' 


(7b) 


Here (^) denotes averaging over fast oscillations. Equa¬ 
tion (|^) describes back reaction of the generated waves 
upon the non-perturbed distribution function fao- In or¬ 
der to calculate the quasilinear term QL it is enough to 
substitute the solution of Eq. (7b) into Eq. ([fa|). The 
solution of Eq. (7b) in the Eourier presentation. 


fai{r, p, t) = J fakip) exp(ikr - iuji^t) dk, (8) 


writes: 
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/ak(p) = - 


Qa 

rUaC . 


dt' exp[ikp — 


Ek ( 1-^ — 

7 ‘^k/ 


(Ek• pO kc 


Wk 


dfaO 

9p' 


( 9 ) 


r 


where Wk = w(k) + zFk; p = r' — r and t = t' — t. 
Then, perturbed electric field vector Ei is substituted 
for perturbed magnetic field Bi using Maxwell equation 
Bk = (c/ci;k)(k X Ek) for the Fourier transforms: 


Ei(r,t) 

Bi(r,t) 


1 

1 


Ek exp(zkr — dk; 


Bk exp(ikr — iuid) dk. 


( 10 ) 


In order to find /Qk(p), we are using standard method 
of integration along non-perturbed trajectories. Follow¬ 
ing this method, in Eq. (^), (r', p') are phase coordinates 
of the particle (along the non-perturbed trajectory) at 
the instant of time t'; they are calculated from the rela¬ 
tivistic equations of motion of a single particle: 


dr' 

dt' 



( 11 a) 


dp' 

dt' 


da 

niaC 


p' xBo{r',t'] 
1 


( 11 b) 


Using cylindrical coordinate system, the non-perturbed 
magnetic field is modeled as Bq = Bo(0, 0, B^p). It is 
supposed that non-perturbed external electric field is ab¬ 
sent in the plasma, Eq = 0. The drift velocity (see Eq. |^) 
is directed along x axis. The particles of opposite charge 
da drift in opposite directions. Hereafter the positive di¬ 
rection of X axis is the direction of the positive charge 
drift. 

In such a geometry, the relativistic equation of motion 
(lit) is rewritten as the following set of equations: 


dt' 7 ’ 


dp'r , dp' 
-Pv 


^-P^-JTT + —Px = 0 


dt 


dt' 

dt' 


1 
/ d^p' 


( 12 ) 


Ignoring the terms proportional to the small parameter 
{tl/RbY (where tl — cjuiBa is the radius of Larmor 


J 


circle) and assuming that p j ^ {p^ +P^)) we can obtain 
the solutions of set of Eqs. (|l2| ) as was done in & 

p'x = Pdc + Pr sin(a;B,:, • r) -t {px - pd„) cos(d;B„ • r), 
p'^ = Pr cos(a;B„ • t) - {p^ - pdJsm{ojBr. ' 't), 

Pip = Pip- (13) 

Here uib^ = ^b^, H\ the components of dimensionless mo¬ 
mentum (pa:, Pn Pp) are the values of (p(,, p)., p(,) at the 
instant of time t' = t, and pd^ = (ud^/c)^. Here we have 
the following integrals of motion: 7 , p(. — loBc'’''! c, r'p'^. 
Let us notice that the particle distribution function fao 
should only depend on the integrals of motion. 

It is evident that the solutions of equation of motion 
( |l^ ) differ from those in the homogeneous magnetic field 
only by the drift term pd^ ■ Using cylindrical coordinates 
in momentum space (py, p±, 9) (subscripts ‘||’ and ‘_L’ 
denote parallel and perpendicular directions to the mag¬ 
netic field Bo, respectively), we can write 


p_L cos 9 = Px-Pda., 

P±sill9 = Pr- 

Therefore Eqs. (|^) are reduced to the following form: 


p'x = Pda.+P±cos{9-UBa.- t), 

p'r = P± sin( 6 * - ujBa. ■ t), 

P'p = P\\- (14) 


Substituting Eqs. & into Eq. 
following expressions for p = r' — r: 


( 11 a) we obtain the 


Px 


Pr 


Pp 


CPdq, 

7 

c 

‘x’Ba. 


Q 

■T -p_L[sin(0 — lvBo, • t) 

UJBa 

p± [cos( 6 * — CoBa.' t) — cos 9], 


^r. 

7 


sin 9 ], 


(15) 


Substituting expressions ( |^ and ( 0 ) for the compo¬ 
nents of p and p' in Eq. (^, we finally obtain Fourier 
transform of oscillating distribution function fauip)- 


/ak(p) = - 


da \ i exp [j 6 sin (0 — c)] 


rriaC 

dfao 

. ^P\\ 


E||(k) + ^^E,(k) 
Wk 7 


-F 


df. 


aO 


dp± 


P\\ OJk 7 P|| 


,( 16 ) 
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where Awak = Wk —fc||'!^|| —kxUd^- Calculating expression 
( 1 ) from Eq. (|^), we used the following presentation of 
the exponential function: 


OO 

exp[zkp — ito-kT] = exp[z 6 sin(d — c)] £n Jn{b)- 

n— — oo 

(17) 


assumptions, the following relations between components 
of electric field and wave vector can be written as 

E\ = 0, Elkx = -ElK (21) 

for purely electromagnetic (E _L k) t waves and 

El^kr = El^kx, = -Ei^kx (22) 


Here = exp[zri(c — 9) — jr(AwQk ~ Jn{b) (n = 

0; ±1; ±2...) is the Bessel function of integer order |M; 


kj_c 

b = -p_L 


(18) 


and c is defined as follows: 


kx = A:^cosc, 

kr = fc^sinc- (19) 

In derivation of Eq. (|l^), we have assumed that 5 <C 1 
and kept only the terms with n = 0 in Eq. (0)- This 
approximation leaves only the terms describing contribu¬ 
tion of Cherenkov-drift resonance. Let us mention that 
{dfao/d9) = 0 since the distribution function possesses 
an axial symmetry. Therefore, corresponding terms do 
not contribute into Eq. (p^. 


IV. QUASI-LINEAR DIFFUSION 


for It waves, which are almost electromagnetic if we as¬ 
sume kx » kr- Using Eqs. ( ^l| ) and (^ we can write 
the diffusion coefficients for t waves and It waves as 


„ = 


n* — 




2 2 

Pd. 


/t(k)dk 


(23a) 


— ^ (1 - —') -^‘(k) (23b) 


/_oo wk r 


Wk 


= Dl ||=0 


for transverse t waves and 


(23c) 



+ ^— Iitik)dk (24a) 



Iitik)dk (24b) 
(24c) 


The next step is to study alteration of slowly varying 
part of distribution function /ao(p,MO due to develop¬ 
ment of Cherenkov-drift instability. Generally, alteration 
is described by diffusion coefficients involved in the quasi- 
linear term. The coefficients show the rate of particle dif¬ 
fusion in momentum space along, as well as across, the 
magnetic field. Substituting Eq. (|^) into Eq. (|^) and 
using Maxwell equation for Fourier transforms, we obtain 
QL term in the following form: 


QL = - 

- A. 

~ 9p\\ 

1 


qa 


rriaC \ dp 

dfaO 


£- I El 


p X B] 


D 


d 


P± dp± 


dpn 

p1d±ii 


d 

^ 9P\\ 

dfaO 

9 p\\ 


7 

P±D\\± 
1 


fal 

dfaO 


P± dp_L 


dp± J 

d f „ dfao \ 


P±D±± 


dp_L J 
( 20 ) 


for longitudinal-transverse It waves. Here we use the fol¬ 
lowing definition: 


It,it(k)=i 



y Eidt{-k)Ei’^\k) 
rriac) V (27rf 


Note that expressions ( p^ ) and ( p^ ) are obtained af¬ 
ter averaging Eq. (^) over the angle 9. This procedure 
nullifies the diffusion coefficients and ; we also 

ignored = 0 taking into account that the terms 


D 


dfao \ 

' 9 pil J 


d 

dp\\ 

are larger then the term 


and 


I _^ 

P± dp± 





d 

dp\\ 


Pi.D\\a 


dfM \ 

dp± J 


Here i4||||, 71||±, and are the diffusion coef¬ 

ficients. Below we calculate particular expressions for 
diffusion coefficients corresponding to t and It waves in 
the case of Cherenkov-drift instability. Let us notice that 
Eq. (p^) has rather general meaning and can as well be 
used for other type of instabilities. However, diffusion 
coefficients will differ for different instabilities. 

It is convenient to consider quasilinear development 
of t and It waves separately. Using simple geometrical 


by a factor (pw/pd.)"^- In Eqs. ( p^ ) and ( p4| ) Ex{—k) 
identifies the complex conjugate to the x component of 
electric field vector Ek, and V = dr. 

The multiplier 

_^ Lk 

^AcJak (^w)^+rk 

7r(5(Aa;), if (Aa;)^ » T^; 
l/Tk, if (Ao .)2 «Tl 
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presenting in the expressions of diffusion coefficients d^ ) 
and ( 1 ^ ) depends on the type of instability: the up¬ 
per case corresponds to the kinetic approximation and 
the lower case to hydrodynamic one. Consequently they 
should be used for instabilities in kinetic and hydrody¬ 
namic approximations respectively. 

It is worth to note that the drift velocity (0) depends 
on the Lorentz factors of the particles, 7 . Hence, ther¬ 
mal spread in Lorentz factors of resonant particles (yr = 
I7—7o|) results in scatter of drift velocities, increasing the 
resonant width of instability, Atu = uj — k\\v\\ —kxUd- It al¬ 
lows to consider the kinetic approximation of Cherenkov- 
drift instability. However, the resonant width of usual 
Cherenkov instability is smaller than the corresponding 
width of Cherenkov-drift instability. Therefore, in the 
presence of narrow relativistic beam (with low value of 
7 t) the kinetic approximation for usual Cherenkov in¬ 
stability is not valid. The growth rate of the insta¬ 
bility, Tk, is small for hydrodynamic approximation as 
well. Therefore usual Cherenkov instability, as opposed 
to Cherenkov-drift instability, can not develop in rela¬ 
tivistic magnetized pair plasma p7| . 

Particle diffusion in momentum space appears both in 
parallel and perpendicular directions with respect to the 
magnetic field Bq. Diffusion process causes alteration of 
particle distribution function until the quasilinear relax¬ 
ation of instability is saturated {dfo/dfit = 0 , where /o 
is distribution function of resonant particles). In order 
to investigate the stage of saturation of quasilinear re¬ 
laxation, it is worth to rewrite Eq. in the following 
form: 


^ = — ( n 
dfit 9p|| V 


1 d 

— o—P-L 
P± Op_L 



(26) 


For estimation of diffusion coefficients D±± and D||||, we 
suppose that kx/k\\ « Pd/j and fc||C w Wk and use the 
following approximations in Eqs. (^,p^: 



Then we rewrite Eq. (p^ in the form of quasilinear inte¬ 


gral of Cherenkov-drift instability: 


dfo 

dpt 


d At^it 
dp\\ [ TT 


1 d 

— ^ 

P± dpj_ 


TT V C / Op± 


(.28) 


Here cuq is the resonant frequency of t and It waves excited 
on Cherenkov-drift instability (see e.g. HI): 

^0 - 3/2^ (7o7p)^^^; (29) 

7t 

At^it is a numerical coefficient (At = 2 and Alt = 8 ); 

denotes the ratio of wave energy and kinetic 

energy of plasma particles (Wp = mc^npjp) for t and It 
waves respectively: 




1^ 


w, 


t.lt 


Wp 

roo TPi-lt 


-dk = 




— 00 
t.lt ( 


mc^npjp V (27r):^ 


^t./i(_k)L;‘.'*(k) 

Stt 


dk; (30) 


the multiplier (i/AuJak) is replaced by l/Tk (p5[), where 
Tk is determined from Eqs. (0) and i) assuming 

{kr/k±) « {kx/k±) w 1. 

Stationary state (dfo/dpt = 0) is reached in the fol¬ 
lowing two cases: 

Pd P±- In this case the first term in Eq. (|^), 
containing alteration of distribution function over 
parallel momenta (dfo/dp\\), is significant. Pro¬ 
cess of quasilinear relaxation will be saturated when 
plateau is formed on parallel distribution function 
of resonant particles (i9/o/9p|| = 0 ); 

t>) Pd ^ P±- In this case the second term of Eq. (|^), 
containing alteration of distribution function over 
perpendicular momenta (dfo/dp±), appears signif¬ 
icant. During the process of quasilinear relaxation, 
energy transfers from parallel to perpendicular mo¬ 
tion of the particles, hence increasing p±. Relax¬ 
ation will be saturated when p± Pd- In this case 
the right hand side terms of Eq. (|^) cancel each 
other. Indeed, they are of the same order but with 
opposite signs (for the beam particle distribution 
function: dfo/dp± < 0 ). 


V. CONCLUSION 

In this paper we discuss the development of Cherenkov- 
drift instability in relativistic magnetized pair plasma. 
We are taking into account particle drift motion across 
the plane of SOME lines, which is significant for the parti¬ 
cles of relativistic beam penetrating the bulk pair plasma. 
We studied quasilinear stage of the instability - quasilin¬ 
ear interaction of excited waves over plasma particles and 















corresponding redistribution of particle momenta. As a 
result, diffusion of particle momenta takes place along, 
as well as across, the magnetic field lines. 

The linear stage of Cherenkov-drift instability devel¬ 
ops similarly to the usual Cherenkov instability. The 
mechanism of wave excitation is based on the well known 
Cherenkov wave-particle interaction. Same as in the case 
of Cherenkov instability, presence of high energy par¬ 
ticle beam (with positive slope on the shape of distri¬ 
bution function) is necessary condition for developing 
of Cherenkov-drift instability. However, in the case of 
Cherenkov-drift instability, the generation of oscillations 
with E T Bo is possible only due to particle drift mo¬ 
tion across Bq (contrary to the case of usual Cherenkov 
instability generating only longitudinal oscillations with 
E II k II Bq). Cherenkov-drift instability generates both 
t (with E±) and It (with E± and E\\) waves; electro¬ 
magnetic oscillations with E^^ are generated by particle 
longitudinal motion with velocities U||, as in the case of 
usual Cherenkov instability. 

Back reaction of excited waves over resonant particles 
should suppress the reason of wave excitation: in the 
case of Cherenkov-drift instability, the process causes, at 
the same time, formation of plateau on the distribution 
function of parallel momenta (similar to the quasilinear 
case of usual Cherenkov instability) and energy transfer 
from parallel motion of particles to their motion across 
the magnetic field. 

The later process, inhibiting anisotropy in momentum 
space, is similar to the quasilinear relaxation of cyclotron 
instability. The reason for development of cyclotron in¬ 
stability - anisotropy in momentum space {p± 'C py) - is 
suppressed by particle diffusion over perpendicular mo¬ 


menta. Perpendicular diffusion is described by nonzero 
diffusion coefficients D±± and D^. As a result, the 
energy of parallel motion of beam particles is transfer¬ 
ring into the perpendicular energy until p± ^ py. As for 
Cherenkov-drift instability, relaxation is saturated since 
Pj_ ~ Pd and the rate of alteration of distribution func¬ 
tion over perpendicular momenta is descr ibed by no nzer o 
perpendicular diffusion coefficients D±± ( ^3b| ) and ( |24tj ). 

It is worth to note that, in the Cherenkov-drift in¬ 
stability, the plateau on the parallel distribution func¬ 
tion forms not only due to parallel diffusion of particles 
(which transfers energy from high velocity resonant par¬ 
ticles to those with low velocities as for usual Cherenkov 
instability), but also because of perpendicular diffusion 
which transfers energy from parallel to perpendicular mo¬ 
tion of particles. This scenario works if the other factors 
which can balance quasilinear diffusion, are not taken 
into account. Such factors could be, on the one hand, 
the radiation reaction force (acting on synchrotron emit¬ 
ting particle, spiraling in strong magnetic field) and, on 
the other hand, the force arising due to particle motion 
in weekly inhomogeneous field. We plan to include these 
factors into consideration in the future works. 
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